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Abstract 

Due to the mathematical complexity of the Boltzmann equation for inelastic hard 
spheres, a kinetic model has recently been proposed whereby the collision rate 
(which is proportional to the relative velocity for hard spheres) is replaced by an 
average velocity-independent value. The resulting inelastic Maxwell model has re- 
ceived a large amount of recent interest, especially in connection with the high 
energy tail of homogeneous states. In this paper the transport coefficients of inelas- 
tic Maxwell models in d dimensions are derived by means of the Chapman-Enskog 
method for unforced systems as well as for systems driven by a Gaussian thermostat 
and by a white noise thermostat. Comparison with known transport coefficients of 
inelastic hard spheres shows that their dependence on inelasticity is captured by 
the inelastic Maxwell models only in a mild qualitative way. Paradoxically, a much 
simpler BGK-like model kinetic equation is closer to the results for inelastic hard 
spheres than the inelastic Maxwell model. 
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1 Introduction 



As is known, the Boltzmann equation for elastic hard spheres is in general 
very complicated to deal with, so that explicit results are usually restricted 
to small deviations from equilibrium [1]. In order to explore a wider range of 
situations, the direct simulation Monte Carlo (DSMC) method [2,3] can be 
used as an efficient tool to solve numerically the Boltzmann equation. From a 
theoretically oriented point of view, another fruitful route consists of replacing 
the detailed Boltzmann collision operator by a simpler collision model, e.g. the 
BGK model [4], that otherwise retains the most relevant physical features of 
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the true collision operator [5]. Several exact solutions of the nonlinear BGK 
model kinetic equation [6] have proven to agree rather well with DSMC results 
for the Fourier flow [7], the uniform shear flow [8], the Couette flow [9,10] and 
the Poiseuille flow [11]. In a third approach, the mathematical structure of 
the Boltzmann collision operator is retained, but the particles are assumed 
to interact via the repulsive Maxwell potential (inversely proportional to the 
fourth power of the distance) [12]. For this interaction model, the collision 
rate is independent of the relative velocity of the colliding pair and this al- 
lows for a number of nice mathematical properties of the collision operator 
[13,14,15]. Since many interesting transport properties (both linear and non- 
linear), once properly nondimensionalized, are only weakly dependent on the 
interaction potential, exact results derived from the Boltzmann equation for 
elastic Maxwell molecules [16] are often useful for elastic hard spheres [8] and 
even Lennard- Jones particles [17,18]. Quoting Ernst and Brito [19,20], one can 
say that "Maxwell molecules are for kinetic theory what harmonic oscillators 
are for quantum mechanics and dumb-bells for polymer physics." 

The prototype model for the description of granular media in the regime of 
rapid flow consists of an assembly of (smooth) inelastic hard spheres (IHS) 
with a constant coefficient of normal restitution a [21]. In the low density limit 
spatial correlations can be neglected. If, in addition, the pre-collision velocities 
of two particles at contact are assumed to be uncorrelated (molecular chaos 
assumption), the velocity distribution function obeys the Boltzmann equation, 
modified to account for the inelasticity of collisions [22,23]. 

Needless to say, all the intricacies of the Boltzmann equation for elastic hard 
spheres are inherited and increased further by the Boltzmann equation for 
IHS, the latter introducing the coefficient of normal restitution < a < 1 
in the collision rule. Therefore, it is not surprising that the three alternative 
approaches mentioned above for elastic collisions, namely the DSMC method, 
the model kinetic equation and the Maxwell model, have been extended to 
the case of inelastic collisions as well. As for the DSMC method, its extension 
from the original formulation for clastic collisions (a = 1) [2,3] to a < 1 is 
straightforward [24]. The generalization of the familiar BGK model kinetic 
equation is less evident, but a physically meaningful proposal has recently 
been made [25]. This model has proven to yield results in good agreement with 
DSMC data for the simple shear flow problem [26,27] and for the nonlinear 
Couette flow [28]. Following the third route, the collision rate of a colliding 
pair of inelastic spheres, which is proportional to the magnitude of the relative 
velocity, is replaced by an average constant collision rate [29]. The resulting 
collision operator shares some of the mathematical properties of that of elastic 
Maxwell molecules and so this model is referred to as pseudo-Maxwellian 
model [29] or, as will be done here, inelastic Maxwell model (IMM) [19]. The 
Boltzmann equation for IMM has received a large amount of interest in the 
last few months [19,20,29,30,31,32,33,34,35,36,37,38,39,40]. 
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The IMM is worth studying by itself as a toy model to exemplify the non- 
trivial influence of the inelastic character of collisions on the physical prop- 
erties of the system. On the other hand, its practical usefulness is strongly 
tied to its capability of mimicking the relevant behaviour of IHS. A common 
property of the Boltzmann equation for IHS and IMM is that both have homo- 
geneous solutions exhibiting high energy tails overpopulated with respect to 
the Maxwcll-Boltzmann distribution. But this general qualitative agreement 
fails at a deeper level. More specifically, in the homogeneous cooling state the 
(reduced) velocity distribution function /(c), where c is the velocity relative 
to the thermal velocity, decays asymptotically as In / ~ — c in the case of IHS 
[41,42] and as ln/~ -Inc in the case of IMM [19,20,33,35]. Analogously, in 
the steady homogeneous state driven by a white noise forcing, the asymptotic 
behaviour is In / -c^/^ in the case of IHS [42], while In / ~ -c for IMM [37]. 
Of course, these discrepancies in the limit of large velocities do not preclude 
that the the IMM may characterize well the "bulk" of the velocity distribu- 
tion of IHS, as measured by low degree moments. In a nonequilibrium gas, the 
physically most relevant moments (apart from the local number density n, flow 
velocity u and granular temperature T) are those associated with the fluxes 
of momentum and energy. If the gradients of n, u and T are weak enough, 
the fluxes are linear combinations of the gradients, thus defining the transport 
coefficients (e.g. the shear viscosity and the thermal conductivity) as nonlin- 
ear functions of the coefficient of restitution. These transport coefficients have 
been obtained by an extension of the Chapman-Enskog method [1] from the 
Boltzmann equation for IHS [43,44,45,46]. To the best of my knowledge, they 
have not been derived for IMM. The major aim of this paper is to carry out 
such a derivation and perform a detailed comparison between the transport 
coefficients for IMM and IHS. 

The plan of the paper is as follows. In Sec. 2 the Boltzmann equation for 
IMM in d dimensions is introduced. The model includes an average collision 
frequency uj that can be freely fitted to optimize the agreement with IHS. In 
the absence of any external forcing the energy balance equation contains a 
sink term due to the collisional energy dissipation. This term is represented 
by the cooling rate that is proportional to the collision frequency ui and to 
the inelasticity parameter 1 — The sink term can be compensated for by an 
opposite source term representing some sort of external driving. For concrete- 
ness, two types of driving are considered: a deterministic force proportional 
to the (peculiar) velocity (Gaussian thermostat) and a stochastic "kicking" 
force (white noise thermostat). The corresponding homogeneous solutions are 
analyzed in Sec. 3, where the fourth cumulant (or kurtosis) of the velocity 
distribution is exactly obtained, being insensitive to the choice of uj. Compari- 
son with the fourth cumulant of IHS shows significant deviations, especially in 
the case of the homogeneous cooling state (which is equivalent to the homoge- 
neous steady state driven by the Gaussian thermostat). The Ghapman-Enskog 
method is apphed in Sec. 4 to get the transport coefficients of IMM in the un- 
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driven well as in the presence of the Gaussian thermostat and the 

white noise thermostat. For undriven systems with c? < 3, it is found that the 
thermal conductivity diverges at a = — d)/3d and becomes negative for 
q; < (4 — d)/3d, irrespective of the choice of uj. A critical comparison with the 
transport coefficients of IHS is carried out in Sec. 5. The free parameter lo is 
fixed by the criterion that the cooling rate of IMM be the same as that of 
IHS (in the local equilibrium approximation). The comparison shows that the 
IMM retains only the basic qualitative features of the a-dependence of the 
IHS transport coefficients. The best agreement takes place in the case of the 
white noise thermostat, where the influence of the inelasticity on the values 
of the transport coefficients is rather weak. Quite surprisingly, the transport 
coefficients predicted by the much simpler BGK-like model [25] are in general 
much closer to the IHS ones than those obtained from the IMM. The paper 
ends in Sec. 6 with some concluding remarks. 



2 Inelastic Maxwell models 

The Boltzmann equation for inelastic Maxwell models (IMM) [20,29,31] can 
be obtained from the Boltzmann equation for inelastic hard spheres (IHS) by 
replacing the term \g-d-\ in the coUision rate (where g = Vi — V2 is the relative 
velocity of the colliding pair and a is the unit vector directed along the centres 
of the two colliding spheres) by an average value proportional to the thermal 
velocity Vq = ^2T/'m (where T is the granular temperature and m is the mass 
of a particle) . The resulting Boltzmann equation is [20] 

(a, + V. . V + ^)/(r, V, t) ^ ^ / da / dv, (a-S- - l) 

x/(r,Vi;t)/(r,V2;t) 
^J[r,^r,■t\f], (1) 

where n is the number density, uj{a) oc nT^^"^ is an effective collision frequency, 
Vtd = 27r°'/^/r((i/2) is the total solid angle in d dimensions, a < 1 is the coef- 
ficient of normal restitution, and b is the operator transforming pre-coUision 
velocities into post-coUision ones: 

-> 1 + a . ^ 

&vi,2 = Vi,2T (g-<r)<r. (2) 

Equations (1) and (2) represent the simplest version of the model, since the 
collision rate is assumed to be independent of the relative orientation between 
the unit vectors g and a [20]. In a more realistic version, the collision rate has 
the same dependence on the scalar product g • o" as in the case of hard spheres. 
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The corresponding Boltzmann equation can be proved to be equivalent to Eq. 
(1), except that the operator b must be replaced by [29,31] 

&vi,2 = 2 (^1 + ^2) ± ± ~^9(r. (3) 

Both versions of the model yield similar results in issues as dehcate as the 
high energy tails [20,37]. For the sake of simplicity, henceforth I will restrict 
myself to the version of the model corresponding to the conventional collision 
rule (2). 

The collision frequency u;{a) is a free parameter of the model. Its detailed 

a-dependence can be determined by optimizing the agreement between the 
results derived from Eq. (1) and those derived from the original Boltzmann 
equation for IHS. Of course, the choice of a;(Q;) is not unique and may depend 
on the property of interest. 

In Eq. (1), is an operator representing the action of a possible external 
driving. This operator is assumed to preserve the local number and momentum 
densities, i.e., 

I dv Tf{v, y;t)^l dv v^/(r, v; t) = 0. (4) 

On the other hand, in general the external driving gives rise to an energy 
source term 



where 



dn{r, t) 

defines the granular temperature and V = v — u is the pecuhar velocity. 



being the flow velocity. A possible external driving corresponds to a determin- 
istic nonconservative force of the form |m7V [37,47] . It has been widely used 
to generate nonequilibrium steady states in the context of molecular fluids and 
can be justified by Gauss's principle of least constraints [48,49]. The operator 
describing this force is 

^/(r,v;i) = ^7(r,t)|^-[V/(r,v;i)]. (8) 

The most commonly used type of driving for inelastic particles consists of a 
stochastic force in the form of Gaussian white noise [30,31,33,42,50,51,52,53,54,55]. 
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Its associated operator is 



The macroscopic balance equations for the local densities of mass, momentum 
and energy follow directly from Eq. (1) by taking velocity moments: 

+ ^^V • u = 0, (10) 



Au+— V-P = 0, (11) 
mn 

In these equations, A = + u • V is the material time derivative, 

P(r, t)^mj dv VV/(r, v; t) (13) 

is the pressure tensor, 

q(r,t) = ^_/ dvl^V/(r,v;t) (14) 

is the heat flux, and 

((r, t) = --— 4?- r / dv V^J[r, v; t\f] (15) 

' ^ dn{r,t)T{r,t) J L ' ' i-'J v ^ 

is the coohng rate. The energy balance equation (12) shows that the existence 
of a driving with the choice 7 = C compensates for the cooling effect due to the 
inelasticity of collisions. In that case, the driving plays the role of a thermostat 
that makes the macroscopic balance equations (10)-(12) look like those of 
a conventional fluid of elastic particles. On the other hand, the transport 
coefficients entering in the constitutive equations are in general different from 
those of a gas of elastic particles and also depend on the type of thermostat 
used. In what follows, I will assume that either 7 = (undriven system, 
JF = 0) or 7 = ^ in Eq. (8) (Gaussian thermostat) and Eq. (9) (white noise 
thermostat). 

The balance equations (10)-(12) arc generally valid, regardless of the details of 
the model for inelastic collisions. However, the influence of the collision model 
appears through the a-dependence of the cooling rate. In the case of IMM, 
one can easily prove (cf. Appendix A) the following relationship between the 
collision frequency cu and the cooling rate (: 

da) = ^^(«)- (16) 
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3 Homogeneous states 



Before solving the inhomogeneous equation (1) by the Chapman-Enskog method, 
it is necessary to analyze the homogeneous solutions, especially their devia- 
tions with respect to the Maxwell-Boltzmann distribution as characterized by 
the fourth cumulant 



where 



{v')^- [ dv^;V(v). 

Th J 



3.1 Homogeneous cooling state. Gaussian thermostat 



In the absence of any external driving (J^ = 0, 7 = 0), Eq. (12) for homo- 
geneous states reduces to dfT — —C,T. It is convenient to scale the velocities 
with respect to the thermal velocity VQ{t) — ^j2T{t)/m and define the scaled 
quantities 

/(c,r)=n-i^;o'W/(v,t), c = Y/vo{t), dr = a;dt. (19) 



Thus, Eq. (1) reduces to 



= Jlci\J], (20) 

where use has been made of Eq. (16). It is interesting to remark that Eq. (20) 
coincides with the Boltzmann equation corresponding to a homogeneous steady 
state driven by the operator (8) with 7 = C (Gaussian thermostat). In other 
words, the application of the Gaussian thermostat to a homogeneous system 
is equivalent to a rescaling of the velocities in the freely cooling case. 

The so-called homogeneous cooling state (HCS) is characterized by a similar- 
ity solution in which all the time dependence of / occurs through the scaled 
velocity c, so that it corresponds to a stationary solution of Eq. (20). Such a 
solution exhibits an overpopulated high energy tail of the form /(c) ~ ^-d-a 
[19,20,33,35], where, in general, the exponent a depends on the coefficient of 
restitution. On the other hand, the "bulk" properties of the velocity distribu- 
tion function are associated with low degree moments. Of course, by normal- 
ization (c^) = d/2. Thus, the first non-trivial moment is (c^). Let us multiply 
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both sides of Eq. (20) by cf and integrate over Ci to get 

dr{c^) - (1 - --l^s = j dci c\J[c,\f\. (21) 

This hierarchy of moment equations can be solved sequentially and has been 
analyzed in detail by Ernst and Brito [20]. For s = 2 we get the identity 
li2 = (1— tt^)/4, which is not but Eq. (16) in dimensionless form. The coUisional 
moment 114 is evaluated in Appendix A with the result 



1^4 



[a{c^) [3a^{l -a)- a{17 + Ad) + 3(3 + Ad) 



l + a 
32d{d + 2) 
-{l + a)d{d + 2){4:d-l-6a + 3a'^)] 



(22) 



Inserting this into Eq. (21) with s = 4 we get the evolution equation 



Ad - 7 + 3a{2 



a 



8d{d + 2) 



d{d + 2) M 



3a{2 



a 



' 4 4d - 7 + 3q;(2 - a) 

(23) 

In the one-dimensional case {d = 1), the fourth moment diverges as (c^) ~ 
exp[r(l — «^)^/8)]. This is consistent with the fact that in this case the exact 
stationary solution to Eq. (20) is /(c) = {2?'^ /ti){1 + 2c^)-^ [19,20,33,35]. On 
the other hand, for d > 2 the moment (c^) relaxes to a well defined stationary 
value (HCS value) whose associated cumulant is 



02 



d{d+2) 



ic') - 1 



a 



Ad - 7 + 3a{2 - a) 



(24) 



This expression is exact for the IMM given by Eq. (1). In contrast, the cumu- 
lant 02 for IHS in the HCS is not known exactly. Nevertheless, an excellent 
estimate is [42] 

^iHS ^ 16(l-a)(l-2a^) 
^ 9 + 2Ad-a{Al-8d) + 30a^{l-a)' ^ ' 

Figure 1 compares the result (24) for IMM with the estimate (25) for IHS. It 
can be observed that the HCS of IMM deviates from the Maxwell-Boltzmann 
distribution (which corresponds to 02 = 0) much more than the HCS of IHS 
[56]. This is consistent with the fact that the former models have a stronger 
overpopulated high energy tail [19,20,33], /(c) ~ c-'^-", than the latter [42], 
/(c) ~ e— . 
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Fig. 1. Plot of the fourth cumulant 02 in the HCS for d = 3 (sohd hnes) and d = 2 
(dashed hnes). 

3.2 White noise thermostat 



Now we assume that the system is heated with the white noise thermostat 
represented by the operator (9) with 7 = C- Using again the scaled quantities 
(19), Eq. (1) becomes 



8d Uci 



/(Ci) = J[Ci|/]. 



(26) 



Taking moments we get 



dAc')-^—^s{s + d-2){c^ 



.s-2\ 



8d 



-pi's 



(27) 



In particular, setting s = 4, 



9,(c^) = -(! + «) 



9 + 12d- a{17 + Ad) + 3a'^{l - a) 
8d(d + 2) 



X 



4 d{d + 2) (3-a)(4rf + 5-3a(2 + a)) 
' 4 9 + 12d - a{17 + Ad) + 3q;2(1 - a) 



(28) 



9 




0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



a 

Fig. 2. Plot of the fourth cumulant 02 in the case of a system heated by a white 
noise thermostat for d = 3 (soHd Hnes) and d = 2 (dashed hnes) . 

In the case of this thermostat, the moment (c^) relaxes to a steady state value 
for any dimensionality d. The corresponding exact expression for the fourth 
cumulant is 

6(1 - a)2(l + a) 



02 = 



9 + 12d- a{17 + Ad) + 3a^{l - a) ' 



(29) 



The cumulant 02 for IHS in the nonequihbrium steady state driven by a white 
noise thermostat can be estimated to be [42] 



JHS 



16(1-q;)(1-2q;2) 



73 + 56d - 3q;(35 + 8d) + 30a^{l - a) ' 



(30) 



Figure 2 shows that the differences between the values of 02 for IMM and 
IHS are less dramatic than in the HCS. The behaviour observed in Fig 2 
is in agreement with the high energy tails /(c) ~ e""'^ for IMM [37] and 
/(c) ~e-«^'^' for IHS [42]. 

It is interesting to remark that, according to Eqs. (24) and (29), 02 oc (1 —a)^ 
for IMM in the small inelasticity limit [29,31], while 02 oc (1 — a) for IHS in 
the same limit [42]. 
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4 Transport coefficients 



The standard Chapman-Enskog method [1] can be generahzed to inelastic col- 
hsions to obtain the dependence of the Navier-Stokes transport coefficients on 
the coefficient of restitution from the Bohzmann equation [43,45,46] and from 
the Enskog equation [44]. Here the method will be applied to the Boltzmann 
equation (1) for IMM. 

In the Chapman-Enskog method a factor e is assigned to every gradient op- 
erator and the distribution function is represented as a series in this formal 
"uniformity" parameter, 

/ = /(°) + 6/(^) + e'f('^ + ■■■. (31) 

Use of this expansion in the definitions of the fluxes (13) and (14) and the cool- 
ing rate (15) gives the corresponding expansion for these quantities. Finally, 
use of these in the hydrodynamic equations (10)-(12) leads to an identification 
of the time derivatives of the fields as an expansion in the gradients, 

9, = af)+e9f) + e29f) + --- . (32) 

In particular, the macroscopic balance equations to zeroth order become 

af^n^O, 9f)u = 0, 9f = -(C - 7)T. (33) 

Here we have taken into account that in the Boltzmann equation (1) the 
effective collision frequency u oc nT^^^ is assumed to be a functional of / only 
through the density n and granular temperature T. Consequently, u^^^ — u, 
ujW = tj(2) = . . . = and, using Eq. (16), C^°^ = C, C^^^ = C^^^ = ■ ■ ■ = 0. It 
must be noticed that, in the case of IHS, (^^^ ~ and (^"^^ is small [43]. 

To zeroth order in the gradients the kinetic equation (1) reads 

■ V + ^) f'Ky) = J[v|/W], (34) 

where use has been made of the properties 

d?^f'\V) = -iC - 7)T^/(°)(V) = • V/W(V), (35) 

the last equahty following from the fact that the dependence of /(°^ on the 
temperature is of the form /^°^(V) = nv^'^f'^^^iy/vo). In the undriven case 

(JF = 0, 7 = 0), as well as in the case of the Gaussian thermostat (8) with 7 = 
C, Eq. (34) is equivalent to the homogeneous equation analyzed in Subsection 
3.1. Therefore, f^^\\) is given by the stationary solution to Eq. (20), except 
that n — > n(r, t) and T — > T(r, t) are local quantities and v — > V = v — u(r, t). 
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Analogously, in the case of the white noise thermostat (9) with 7 = C Eq. (34) 
is equivalent to the homogeneous equation analyzed in Subsection 3.2. 

Since /(°^ is isotropic, it follows that 

P(°)=pl, q(°) = 0, (36) 

where p — nT is the hydrostatic pressure and I is the d x d unit tensor. 
Therefore, the macroscopic balance equations give 

= -nV • u, Df^u=- — , Df¥=-^V-u, (37) 

mn d 

where D^^^ = dj:^^ + u ■ V. To first order in the gradients Eq. (1) leads to the 
following equation for /^^^: 

{dP +C + J^) /(^)(V) = - {Di'^ + V • V) /(°)(V), (38) 

where £ is the linearized collision operator 

X /<">(Vi)/<»(V2) + /l»>(V2)/W(V,)] . (39) 

The collisional integrals of VV and K^V are evaluated in Appendix A. Prom 
the linearization of Eqs. (A. 7) and (A. 9) we have 



m 



y dVVV£/(^)(V) =i/P(^), (40) 



^ / dVl/V£/W(V) = (^-+ ^C) q^^\ (41) 
where the collision frequency 1/ is 

d{d + 2) ^ ' 

Using (37), the right-hand side of Eq. (38) can be written as 

-(d^^ + V- V)/(°)(V) = A(V)-VlnT + B(V)-Vlnn + C(V) : Vu, (43) 
where 

B ^ -V/(°) - -j^/<»), (45) 
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Now we multiply both sides of Eq. (38) by mViVj and integrate over V. The 
result is 

(af + u)p^^ + ng) = -pA,,,iVkUi, (47) 

where 

njj) = m| dvy,y,.F/(i)(v) (48) 

and 

2 

AjjM = 5jfc5jz + SiiSjk - -^SijSki- (49) 
Of course, Hj-j-* = in the undriven case. For driven systems, 

en —CP-i^P (Gaussian thermostat), 
ng^ = <^ ^ ^ (50) 

I (white noise thermostat). 

The solution to Eq. (47) has the form 

= -TjAijkiVkUi, (51) 
where r] is the shear viscosity. By dimensional analysis, r] oc T^/^. Therefore, 

a(o)p{i)^_C^p(i). (52) 

Consequently, Eq. (47) yields 



(^u — (undriven system), 

(z/ — C)~^ (Gaussian thermostat), (53) 
(white noise thermostat). 



T) = p < 



Except for a possible a-dependence of [cf. Eq. (68) below], the result rj = p/u 
in the case of the white noise thermostat is the same as for elastic particles. 
This allows us to interpret z/~^ as the effective mean free time associated with 
momentum transport. This formal equivalence between the shear viscosity of 
a fluid of clastic particles and that of a granular fluid driven by a white noise 
forcing is due to the fact that the latter forcing, while compensating for the 
inelastic cooling, does not contribute to the rate of change of the stress tensor. 
The Gaussian thermostat, on the other hand, yields a term CP^^-* and therefore 
tends to produce a temporal increase in the magnitude of the stress tensor, 
thus partially cancelling the dissipativc term —uP^^\ As a consequence, the 
steady-state shear viscosity is enhanced, rj = p/{i^ — C)- Finally, in the absence 
of any external driving, the state is unsteady and so the inelastic cooling is 
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responsible for a smaller enhancement of the shear viscosity, f] = p/{v — C/2). 
It is worth noting that a structure similar to that of Eq. (53) is also present in 
the cases of the Boltzmann equation for IHS [cf. Eq. (B.l)] and the BGK-like 
kinetic model [cf. Eq. (C.3)]. 

Let us consider next the heat flux. Multiplying both sides of Eq. (38) by 
^mV^Y and integrating over V we get 

•2 



(af ) + u') qW + QW = -^(1 + 2a,)^VT - ^a.^Vn, (54) 
^ ^ 2 m 2 m 

where 



, d-1 d + 4(d-l) + (8 + d)(l-a) 

u = — ; — 1/ H — C = ^ — -\ -ly 55 

d 2d ^ 4d + 4(l-Q;) ^ ^ 

is an effective collision frequency associated with the thermal conductivity and 

Q(i) = dVy2v^/(^)(V). (56) 
In the undriven case, Q*^^-* = 0. For the types of driving we are considering. 



Q 



(1) 



— |Cq^^^ (Gaussian thermostat), ^^^^ 
(white noise thermostat). 

The heat flux has the structure 

q(i) = -AVT - /xVn, (58) 

where A is the thermal conductivity and // is a transport coefficient with no 
counterpart for elastic particles [43] . Dimensional analysis shows that A oc T^/^ 
and II oc T^/^. Consequently, 



di'^q^'^ = 1{C- 7)AVT + ^(C - 7)/^Vn + AV(C - 



= (C-7) 



3 T 
2AVT +(-// + A 



n 



Vn 



(59) 



where in the last step we have taken into account that C; 7 c< nT^^^. Inserting 
this equation into Eq. (54), we can identify the transport coefficients as 



pd+2. ^ . 
A = — ^^(l + 2a2) < 



m 



{v' — 2C,)~^ (undriven system), 

— \C) (Gaussian thermostat), (60) 
(white noise thermostat). 



/-I 
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A 



nl + 2a2 



{( + a2i'') {y' — (undriven system), 

(Gaussian thermostat), 
(white noise thermostat) . 



02 

02 



(61) 



In Eqs. (60) and (61), the cumulant 02 is given by Eqs. (24) (undriven system 
and Gaussian thermostat) and (29) (white noise thermostat). 

Using Eqs. (16) and (42), it can be seen that the thermal conductivity in the 
undriven case is A oc (a — ao)"''", where ao = (4 — d)/3d. This imphes that the 
coefficients A and fi exhibit an unphysical behaviour for d = 2 and d = 3 since 
they diverge at a = ao and become negative for < a; < ao- This singular 
behaviour is absent in the shear viscosity or in A and /j, for thermostattcd 
states because u > ( ior all a and d. 

In order to have the full a-dependence of the transport coefficients we need 
to fix the free parameter a; (a). This point will be addressed in Section 5. On 
the other hand, the ratios between transport coefficients are independent of 
the criterion to choose u!{a). Let us define a generalized Prandtl number (for 
d^l) 



Rn{a) 



A(a)/Ao' ^^^^ 

where 770 and Aq are the shear viscosity and thermal conductivity, respectively, 
in the elastic limit (a = 1). From Eqs. (53) and (60) we have 



1 + 2a2 



d{d- l+:j<ia) 



-l)[3d+2+((i-2)c 



d[5d+4-(d+8)a] 
4(d-l)(d+l-a) 



Analogously, we can define the ratio 



Thus, 



(undriven system), 
(Gaussian thermostat), 
(white noise thermostat) . 



nfi{a) 



T\{a)' 



(63) 



(64) 



1 + 2a2 



a2 (Gaussian thermostat), 

02 (white noise thermostat). 



(65) 
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5 Compcirison with the transport coefficients of inelastic hard spheres 



The transport coefficients of IHS described by the Boltzmann equation have 
been derived both for undriven [43,44,45] and thermostatted [46] systems in 
the first Sonine approximation. For the sake of completeness, the expressions 
of the transport coefficients of IHS are listed in Appendix B. 

Figures 3 and 4 compare the ratios Rr,{a), Eq. (62), and Rn{oi), Eq. (64), for 
IMM and IHS in the three-dimensional case. It is observed that, in general, 
the results for IMM describe qualitatively the a-dependence of i?^ and for 
IHS. Thus, as the inelasticity increases, the generalized Prandtl number i?^ 
decreases in the absence of external forcing and increases in the case of the 
white noise thermostat. With the Gaussian thermostat, however, the discrep- 
ancies arc important: R^^ increases with the inelasticity for IHS and decreases 
for IMM. As for the ratio i?^, which measures the new transport coefficient 
relative to the thermal conductivity, it rapidly increases in the unforced case, 
while it is very small in the driven cases. At a quantitative level, the IMM 
results exhibit important deviations from the IHS ones, especially in the un- 
driven case, where i?^ for IMM becomes negative when « < = | and R^ 
grows too rapidly. The situation in which the IMM ratios are the closest to 
the IHS ones corresponds to the system heated by a white noise thermostat. 

Of course, the most interesting comparison refers to the three transport coeffi- 
cients themselves, rather than to their ratios. In order to have explicit expres- 
sions for the transport coefficients of IMM, we now need a criterion to fix the 
free parameter uj{a). The most natural choice to optimize the agreement with 
the IHS results is to guarantee that the coohng rate for IMM, Eq. (16), be 
the same as that for IHS. Strictly speaking, the cooling rate for IHS depends 
in general on the details of the nonequilibrium velocity distribution function, 
while the effective collision frequency uo{a) in Eq. (1) is assumed to depend 
on / only through the density and the granular temperature. Otherwise, the 
complete knowledge of uj would require to solve first the Boltzmann equation 
for IHS and then evaluate the cooling rate associated with such a sohition, 
what is impractical. Therefore, here I take for C, the cooling rate of IHS at 
local equilibrium, namely, 

C = ^o^(l-«^'), (66) 

where is given by Eq. (B.4). Making use of Eqs. (16) and (42), this is 
equivalent to 

ZiOj 
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0.0 0.2 0.4 0.6 0.8 1.0 

a 

Fig. 3. Plot of the generalized Prandtl number Ii,j{a) from the Boltzmann equation 
for three-dimensional IMM (solid lines) and IHS (dashed lines) for (a) undriven 
systems, (b) Gaussian thermostat and (c) white noise thermostat. The curves cor- 
responding to IMM in case (b) and to IHS in case (a) are very close each other by 
accident. 



With this choice, Eqs. (53), (60) and (61) become, respectively. 



2d 



V = Vo 



1 + a 



4 [3d + 2 + {d — 2)a] ^ (undriven system), 

I (1 + q;)~^ (Gaussian thermostat), 



(69) 



{d + l-a)-^ 



(white noise thermostat). 



l{d-l) 
1 + a 



(1 + 2a2 



(d - 4 + 3da)~^ 

[5d + 4 — {d + 8)a]~^ (white noise thermostat), 

(70) 



(undriven system), 
(Gaussian thermostat). 
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0.0 0.2 0.4 0.6 

a 



0.8 1.0 



Fig. 4. Plot of the ratio Rn{a) from the Boltzmann equation for three-dimensional 
IMM (solid lines) and IHS (dashed lines) for (a) undriven systems, (b) Gaussian 
thermostat and (c) white noise thermostat. 



At = — Ao 



n (1 + a) 



2 



+ (™<biven system), 

a2 (Gaussian thermostat), 

2((j-l)(l+a) 



5d+4-{d+8)a ^"^ (white noise thermostat), 

(71) 

Figures 5-7 compare the three transport coefficients of IMM with those of 
IHS. For completeness, also the coefficients derived from a simple BGK-like 
model (see Appendix C) arc included. Again, the qualitative behaviour of IHS 
is generally captured by the IMM. We observe that the shear viscosity in- 
creases with the inelasticity, the increase being more (less) important when 
the system is heated with a Gaussian (white noise) thermostat. The thermal 
conductivity increases with the inelasticity in a significant way in the undriven 
case, increases more moderately in the case of the Gaussian thermostat, and 
is almost constant in the case of the white noise thermostat. As for the trans- 
port coefficient /i, it remains small in the thermostatted states, while it rapidly 
increases in the undriven state. All these trends are, however, strongly exagger- 
ated by the IMM in the undriven case (where A and // diverge at o; = ccq = \) 
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a 

Fig. 5. Plot of the reduced shear viscosity ??/r/o from the Boltzmann equation for 
three-dimensional IMM (solid lines) and IHS (dashed lines), as well as from the 
BGK-like model of IHS (dotted lines), for (a) undriven systems, (b) Gaussian ther- 
mostat and (c) white noise thermostat. The dashed and dotted lines are practically 
indistinguishable in case (c). 

and, to a lesser extent, in the Gaussian thermostat case, especially for the 
coefficient /x. On the other hand, the simple BGK-likc model [25] summarized 
in Appendix C describes fairly well the a-dependence of the three transport 
coefficients, being closer to the IHS results than the IMM predictions. 



6 Concluding remarks 

In the Boltzmann equation for inelastic Maxwell models (IMM), the collision 
rate of the underlying system of inelastic hard spheres (IHS) is replaced by 
an effective collision rate independent of the relative velocity of the colliding 
particles. Based on the experience in the case of elastic particles, one might 
reasonably think that this model, while making the collision operator mathe- 
matically more tractable, is able to capture the most important properties of 
IHS, at least those relatively insensitive to the domain of velocities much larger 
than the thermal velocity. In particular, one might expect the three transport 
coefficients (shear viscosity, thermal conductivity and the coefficient associ- 
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a 

Fig. 6. Plot of the reduced thermal conductivity A/Aq from the Boltzmann equation 
for three-dimensional IMM (solid lines) and IHS (dashed lines), as well as from 
the BGK-like model of IHS (dotted lines), for (a) undriven systems, (b) Gaussian 
thermostat and (c) white noise thermostat. Note that the curves corresponding to 
IMM and the BGK-like model in case (c) are hardly distinguishable. 

ated with the contribution of a density gradient to the heat flux) of IMM to 
possess a dependence on the coefficient of restitution similar to that of the 
transport coefficients of IHS. The results derived in this paper show, however, 
that this expectation does not hold true, except at a mild qualitative level, 
even if the IMM is made to reproduce the cooling rate of IHS. On the other 
hand, a model kinetic equation based on the well-known BGK model for elastic 
collisions presents a surprisingly good agreement with IHS, once the limitation 
of the BGK model to reproduce the correct Prandtl number of elastic particles 
is conveniently accounted for. A partial explanation of this paradox may lie in 
the fact that in the IMM a rough approximation (collision rate independent 
of the relative velocity) coexists with the remaining complexity of the detailed 
collision process. As a result, only a free parameter (essentially the cooling 
rate) is available to make contact with IHS. In the BGK-Hke model, however, 
the cooling rate and its dominant effect on the velocity distribution function 
is explicitly taken out; what remains of the Boltzmann collision operator is 
modelled by a conventional relaxation-time term reflecting the effects of col- 
lisions not directly associated with the energy dissipation. As a consequence. 
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0.0 0.2 0.4 0.6 0.8 1.0 

a 

Fig. 7. Plot of the reduced transport coefficient ^/{TX^/n) from the Boltzmann 
equation for three-dimensional IMM (solid lines) and IHS (dashed lines), as well 
as from the BGK-like model of IHS (dotted line), for (a) undriven systems, (b) 
Gaussian thermostat and (c) white noise thermostat. Note that the BGK-like model 
predicts /x = in cases (b) and (c). 

in addition to the cooling rate, the BGK-like model incorporates an effective 
(velocity-independent) collision frequency that increases the flexibility of the 
model in spite of its simplicity. 

In conclusion, the IMM is interesting as a mathematical toy model to explore 
how a small degree of coUisional inelasticity may have a strong influence on 
the physical properties of the system. For instance, the velocity moments in 
homogeneous states can be evaluated only approximately in IHS [42], while 
they can be exactly obtained in IMM [29]; the high energy tail is another ex- 
ample where rather detailed information can be obtained from the Boltzmann 
equation for IMM [19,20,33,35,37]. On the other hand, if one is looking for a 
"shortcut" to know some of the properties of IHS, the use of the IMM requires 
a great deal of caution. First, the solution of the Boltzmann equation for IMM 
may still be a formidable task, except in some special cases. Second, the results 
derived from the IMM may not be sufficiently representative of the behaviour 
of IHS, as exemplified here in the case of the transport coefficients. From that 
point of view, it seems preferable to make use of the BGK-hke model proposed 



21 



in Ref. [25]. For instance, this model has an exact solution for the nonlinear 
planar Couette flow (with combined momentum and energy transport) that 
compares quite well with DSMC data for IHS [28]; for this flow, however, the 
Boltzmann equation for IMM cannot be solved in a closed form, even in the 
elastic limit [16]. 
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A CoUisional moments 

Let us consider the general collisional integral of the form 



Now we particularize to -F(v) = VV. Prom the collision rule (2) it follows 
that 




(A.l) 



By following standard steps, I[F] can be written as 



I[F] 



2nQd 



I dvi I dv2 1 da/(vi)/(v2)(6-l)[F(vi)+F(v2)]. (A.2) 



(b-l) (V1V1+V2V2) 



1 + a 



(g-a-) [(1 + a){g ■ a)aa - go" - 6-g] . (A.3) 



2 



To perform the angular integrations we need the results 




(A.4) 




(A.5) 



where [20] 




(A.6) 



and I is the dx d unit tensor. Therefore, 
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/|VVl = -u^ j dv,y"dv.2/(vi)/(v2) 



1 + a 

— — r 
2ma 



(l-a)pl + 2 ^ + |^" (P-pl) 



(A.7) 



The cooling rate is C = — (m/dnTj/fF^], so taking the trace of Eq. (A.7) we 
get Eq. (16). 

Next we particularize to -F(v) = V'^Y. The coUision rule gives 



(h - 1) {v,'y, + viv,) = (^)' (g • ^)'(Vi + V2) • (I + 2aa) 



2 

V^\ - 2V2V2) • a. (A.8) 



After performing the angular integrations one has 



/[^'V] = 5^^^!^ / ^^^/ dV2/(Vi)/(v2) 



X 



-g • (T/i^l + 2V1V1 - Vi\ - 2V2V2)} 

2u;{l + a){d- 1) ^ d + 81-a \ ^ 



md{d + 2) \ d-1 A 

Finally, let us evaluate the coUisional moment where Hs is defined by the 
last identity in Eq. (21), i.e., 

= / d^ / dci J dc2 /(ci)/(c2) (6 - 1) (4 + 4). (A.IO) 

Prom the collision rule one gets [42] 



(b - 1) (4 + 4) = 2(1 + «)'(g • ^)'(G ■ + ^(1 - «')'(g • ^) 

-(l-a^)(g-a)^ (G^ + \g 



4 



-4(l + a)(g-^)(G-^)(g-G), (A.ll) 
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where G = |(ci + C2). Equations (A. 4) and (A. 5) give 



/da(g.a«G.^f = -^[2(g.G)' + ,^G 



\2 I „2/-(2 



J d^(g-a)(G-^) = ^(g-G) 



(A.12) 
(A. 13) 



Therefore, 



//4 = -^^^T- / f ciC2 /(Ci)/(C2) J 



2d 



a 



{ d + 2 



2{g-Gf + g'G' 



- + - (1 - a)g^ (g^ + \9')- 4(g ■ G)^ } (A.14) 



^ 8(ci+2) 
Finally, by taking into account that 



{g-Gr^-i4 + 4-2clcl), 



g'G' = \[4 + 4 + 2clcl-4{c,-c,n 



g^ = 4 + 4 + 244 + 4(ci • C2)' - 4(4 + C2)(Ci ■ C2), 



(A.15) 

(A.16) 
(A.17) 



we get Eq. (22). 



B Transport coefficients of inelastic hard spheres 



The derivation of the Navier-Stokes transport coefficients for IHS can be found 
in Refs. [43,44,45,46]. In the first Sonine approximation, the results are 



p_ 



{v* — \C*^ (undriven system), 
{v* — C*j (Gaussian thermostat), 
(white noise thermostat) , 



(B.l) 



A 



p d + 2 
muQ 2 



(1 + 2ar) 



{ul — 2C,*) ^ (undriven system), 

(^A ~ |C*) ^ (Gaussian thermostat), (B.2) 



(white noise thermostat). 
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T A 
n 1 + 2a^HS 



(C* + ^^^^^a) - fC*) (undriven system), 

(Gaussian thermostat), 



4"^ 



In these equations [45] , 



(white noise thermostat). 

(B.3) 



uq = na (T/m 



^l/2_ 



7rV2((i + 2)' 



IHS 
0^2 I ) 



(B.4) 

(B.5) 
(B.6) 



1 + a 
d 



d — 1 3,, , , 4 + 5(i — 3(4 — d)a ttto 

^ (d + 8) 1 - a) H ^ ^4 

2 16^ ^ 512 ^ 



(B.7) 



The expressions for are given by Eq. (25) for the undriven case and for 
the Gaussian thermostat and by Eq. (30) for the white noise thermostat. 



C BGK-like model 



In its simplest version, the BGK-like model proposed in Ref. [25] reads 



{dt + v-V + T)f{r, v;t)^- (z/BGK - C) [/(r, v; t) - /o(r, v; t)] 

+^c|;-V/(r,v;t), (C.l) 

where i^BCK nT^/^ is an effective coUision frequency and /o is the local equi- 
librium distribution function. The BGK-like model (C.l) can be interpreted 
as indicating that a system of inelastic hard spheres behaves essentially as a 
system of elastic hard spheres with a modified a-depcndcnt rate of collisions 
and subjected to the action of a "friction" force — |mCV, which accounts for 
the energy dissipation in an effective way. Applying the Chapman-Enskog 
method to Eq. (C.l) one obtains Eq. (38) with /(°) = /o and the replacement 

>C/(^^-(-BGK-C-^C^-v)/W. (C.2) 
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It is then straightforward to get the transport coefficients as 

(j^BGK — |C) (undriven system), 
V — P { (i^BGK — C)~^ (Gaussian thermostat), 
i^BGK"''^ (white noise thermostat) , 



(C.3) 



, p d + 2 
A = — < 



m 



fjL = —A < 
n 



(z^BGK - |C) (undriven system), 
^BGK — C)^^ (Gaussian thermostat), 
(i^BGK + |C) (white noise thermostat), 

C (^BGK — C)~^ (undriven system), 
(Gaussian thermostat), 

(white noise thermostat). 



(C.4) 



(C.5) 



Strictly speaking, the coohng rate C is a functional of the velocity distribution 
function through Eq. (15). In the spirit of the kinetic model (C.l), the exact 
cooling rate ( is approximated by its local equilibrium expression 



da) ^ 1^0^(1 - a') 



(C.6) 



where uq oc riT^^"^ is a collision frequency (independent of a) defined by Eq. 
(B.4). It remains to fix the a-dependence of i^BGK- Comparison between Eqs. 
(B.l) and (C.3) suggests the identification ubgk i^oi^^ as the simplest choice. 
Thus, 

(C.7) 



3 / 2 \ 

i^BGK = i^oj^ \ l-a + -dj (1 + a). 



where we have set — > 0, in consistency with the local equilibrium approx- 
imation (C.6) for the coohng rate. Inserting Eqs. (C.6) and (C.7) into Eqs. 
(C.3)-(C.5), one gets 



Ad 
1 + a 



2[A + 3d — {4 — d)a]~^ (undriven system), 

[1 + d + {d — l)oi]~^ (Gaussian thermostat), (C.8) 

(3 + 2d — 3a)~^ (white noise thermostat). 



A (a) = Ao 



4d 
1 + a 



2 [d{l + 3a)]~^ (undriven system), 

[1 + d + {d— l)oc\~^ (Gaussian thermostat), (C.9) 

2 [8 + 5d — (8 + dja]"^ (white noise thermostat). 
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T 



n 



-Ao- 



8{d + 2){l 



a] 



1 + a)(l + 3a)[l + + (rf - l)a] 



1 (undriven system), 

(Gaussian thermostat), 

(white noise thermostat), 

(C.IO) 

In the above equations we have taken into account that rjo/rnXo = 2/{d + 2) 
in the BGK model, while the Boltzmann value is r]o/mXo = 2{d— l)/d{d + 2). 
This discrepancy in the Prandtl number is a known feature of the BGK model 
for elastic collisions and is due to the fact that the model contains a single 
relaxation time. 
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